Spin versus Lattice Polaron: Prediction for Electron-Doped CaMnC>3 
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CaMnOa is a simple bi-partite antiferromagnet(AF) which can be continuously electron-doped up 
to LaMnC>3. Electrons enter the doubly degenerate E g subshell with spins aligned to the S — 3/2 
core of Mn 4+ (TjJ). We take the Hubbard and Hund energies to be effectively infinite. Our model 
Hamiltonian has two E g orbitals per Mn atom, nearest neighbor hopping, nearest neighbor exchange 
coupling of the S — 3/2 cores, and electron-phonon coupling of Mn orbitals to adjacent oxygen atoms. 
We solve this model for light doping. Electrons are confined in local ferromagnetic(FM) regions (spin 
polarons) where there proceeds an interesting competition between spin polarization (spin polarons) 
which enlarges the polaron, and lattice polarization (Jahn- Teller polarons) which makes it smaller. 
A symmetric 7-atom ferromagnetic cluster (Mnf + ) is the stable result, with net spin S~2 relative 
to the undoped AF. The distorted oxygen positions around the electron are predicted. The model 
also predicts a critical doping x ~ 0.045 where the polaronic insulator becomes unstable relative to 
a FM metal. 
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I. INTRODUCTION 

CaMn03 is a bi-partite (G-type) antiferromag- 
netic(AF) insulator jlj with Neel temperature Tat=125K, 
and almost perfect cubic perovskite crystal structure. 
There is not a large literature on this material. It de- 
serves attention as a particularly simple case of an AF 
insulator which can be electron-doped. The (T, a;) phase 
diagram of the Cai-^La^MnOs series has attracted 
attention because of the fascinating interplay of spin or- 
der, orbital order, and metallic versus insulating trans- 
port. "Colossal magnetoresistance" (CMR) occuring at 
concentration x w 0.65 and T w 250 K is the most 
dramatic manifestation [Q. For small x, magnetization 
and conductivity measurements |5||| suggest local ferro- 
magnetic (FM) regions or "spin polarons" in the range 
0.02 < x < 0.06. In this paper we use a model for pure 
CaMnC>3, the x — end member, and predict its behav- 
ior under light doping, x -C 1. We keep T equal to and 
neglect lattice zero-point energy, but all other degrees of 
freedom are allowed and interact in interesting ways. 

Our main question is, what is the ground state of an 
excess electron in CaMnOa? The answer to this question 
will also apply to lightly doped CaMnOa provided the 
doped state is homogeneous. Our model has four param- 
eters. (1) The bandwidth parameter t governs the effec- 
tive Mn Eg electron hopping between Mn 3+ and Mn 4+ 
through the intervening oxygen. The spins of the two Mn 
ions must be parallel as in the usual "double exchange" 
model (?]|| . (2) The magnetic exchange parameter J cou- 
ples spins of first-neighbor Mn ions, due to virtual hop- 
ping of T2 g electrons. (3) The electron-phonon coupling 
constant g describes interactions between Mn E g orbitals 
and the 6 nearest oxygen atoms. (4) Oxygen displace- 
ments u are opposed by the restoring force —Ku, where 
lo = \J K/M is an Einstein frequency assigned to oxygen 
vibrations along the bonds. There are two important di- 
mensionless parameters. Spin polarons || are controlled 
by the ratio (3 = t/JS 2 . Jahn- Teller (JT) lattice po- 



larons |hJ are controlled by the parameter T — g 2 /Kt 
p[ . Balancing these competing effects, we find the most 
favorable local FM spin arrangement, lattice distortion 
and electron wavefunction. As doping increases, we pre- 
dict a transition from polaronic insulator to FM metal. 



II. MODEL HAMILTONIAN 

The Mn 4+ ion in CaMn03 has configuration 3d 3 , i.e., 
the three spin-aligned Ti g states (xy,yz,zx) are filled 
with electrons, while the two spin-aligned E g states 
(-02 = (x 2 — j/ 2 )/\/3, ip3 — 3z 2 — r 2 ) are empty and lie 
above by the crystal field splitting. The empty opposite- 
spin T2 g states are split to even higher energy by the 
Hund term Jh- Light electron doping puts carriers into 
the doubly-degenerate E g level. Hopping of (dd7r)-type 
occurs from Ti g to T2 g , and of (cZcfcr)-type from E g to 
Eg, but no E g to T2 g hopping matrix element exists be- 
cause of the simple cubic structure. Virtual T2 g hopping 
]i~3t (at the cost of Hubbard energy U) gains dereal- 
ization energy if adjacent spins are antiparallcl. This 
gives an S = 3/2 antiferromagnetic Heisenberg Hamil- 
tonian with exchange coupling J = 2(dd-n) 2 jU fl^] ; U 
being large compared with (ddir) j^j, and agrees with 
the experimentally observed magnetic structure of pure 
CaMn03. The ground state has f spins on the A sublat- 
tice (when exp(iQ • V) — 1, where I labels the Mn sites), 
and I spins on the B sublattice (when exp(iQ ■ I) = — 1) 
with Q — (jr, 7r, 7r). 

Ignoring for now the electron-phonon terms, the 
Hamiltonian for an excess electron is Ti = Tit + Hj. The 
first term contains hopping of Mn E g electrons to nearest 
neighbors, 



7U =tJ2[S(l,l± z)c\{T± z)c 3 {l) 
i,± 

+ rotations to x,y directions], (I) 
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6(1, 2) = cos — cos — + sm 



— sin— e ^ v 



(2) 



Here C3(Z) destroys an electron state ^(r — Z) on the Mn 
atom at Z; Z± 5 labels the Mn neighbors above and below 
the one at site Z; and t is the (dda) integral from Slater- 
Koster two-center theory [15[. We use a value f = —0.75 
cV obtained from fitting the band structure of CaMn03 
(see Sec. III). The factor 6(1,2) comes from lo- 
cally rotating the axis of spin quantization for the i th E g 
electron into the direction {0i,4>i) of the i th S = 3/2 
core spin, treating the angles (#i, 4>i) as classical parame- 
ters, and discarding from the Hilbert space the state with 
spin opposite to the core spin (i.e. assuming Jh — > oo 
0,^8|). In this paper, we usually take the spins to be 
perfectly ordered at T — 0, that is, 6 — 61 — 62 equals 
or 7r, corresponding to 6(1,2) equal to 1 or for FM or 
AF oriented neighbors. We will also consider uniformly 
canted states where a relative angle 6 = n — 8q occurs, 
with </>=constant, so that 6 takes the value sin(#o/2). 
These spin orientations are shown in Fig. ^. 




FIG. 1. Schematic spin structures for the antiferromag- 
netic G (AFG), ferromagnetic (FM), and antiferromagnetic 
C (AFC) structures, and interpolating canted structures. 

The rotation of ip3 to the x axis is (—"03 + V / 3'02)/2 = 
3x 2 —r 2 , and to the y axis it is (— ^3— v3~V>2)/2 = 3y 2 —r 2 . 
Using these, we rewrite TL t in the usual orthonormal basis 

(02,^3), 



(3) 



where it is understood that the hopping only operates 
between parallel spin Mn atoms. The hopping matrices 
are 



T x = 



3 
4 

vf 
4 




T z = 




1 

(4) 

The Ttj term is the AF nearest neighbor Ti g exchange: 

Hj= JS(l)-S(T) (5) 

<l,l'> 



The exchange coupling J is can be estimated from Mean 
Field theory (we use a quantum treatment for spin 3/2; 
FM and AFM answers are equal) to be JS 2 = 3.23 meV, 
using the measured Neel temperature Tn = 125 K [fl9|| . 
However, for a given coupling J, a more accurate esti- 
mate from susceptibility expansions reduces the Tq of 
the Heisenberg ferromagnet EQ] 



T, 



c 



Mz-i)[us(s + i)-i] 



T, 



C(MF) 



|z6(6- 



1) 



(6) 



where z = 6 is the number of nearest neighbors. For a 
Heisenberg antiferromagnet |2l],|22j it is estimated that 
Tjv is slightly higher, 



^-14 
T c zS(S 



0.63 



1)' 



(7) 



Making these corrections, the coupling is found to be 
JS 2 = 4.74 meV, giving i/J6 2 =158, with probable un- 
certainty of 10%. The phonon parts of the Hamiltonian 
are given in Sec. V. 



III. UNIFORM SOLUTIONS 

First consider the hypothetical case of a uniform FM 
spin order. Then a doped-in electron could hop without 
paying a Hund penalty, and extended Bloch states would 
form with wavefunctions 

* fe = d 2 4.E eltt Ml) + da-^E ^'^(0, (8) 



N 1 " v^* 



where cfe, ^3 are coefficients. Diagonalizing Eq.(|^), the re- 
sulting energy eigenvalue, taking into consideration that 
t is negative, is 



E(k) 



— cos k x — cos k y — cos k 



=F(cos 2 k x + cos 2 k v + cos 2 k. 
— cos k x cos k y — cos k y cos k z — cos k z cos 



'z 

k x ) ^ ■ 
(9) 



At cosk x — cos k y = cosk z — 1, the energy is minimum, 

E(k = 0) = — 3|t|, the state being doubly degenerate. 
Although the actual spin arrangement is not FM, the re- 
sult |3i| sets a useful scale for the maximal energy gain 
from electron derealization. Also, the dispersion relation 
(^) can be fitted to a FM CaMnC>3 majority-spin band 
structure calculated in density-functional theory Jl6| , de- 
termining the hopping energy t = —0.75 eV. 

For light doping, the magnetic energy Eqn.(^|) of anti- 
ferromagnetic order (—zJNS 2 /2 for classical spins; N 
is the number of Mn ions) is much larger than hop- 
ping energy. As discussed by deGennes p3|, the best 
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uniform solution (for a one-band model) is a compro- 
mise where antiferromagnetic spins cant uniformly to- 
ward FM solutions. Start from the Neel AF structure 
and let all spins tilt towards the x direction with angle 
6/2, as shown in Fig. [I]. This costs magnetic energy 
AEj = zNJS 2 {\ — cos#)/2. A doped-in electron can 
now delocalize, reducing its energy by — i\t\ sin(0/2). For 
small doping x — n/N <C 1, the total hopping energy 
is AE t = -3n|t|sin(0/2). Therefore when sin(0/2) = 
x\t\/4JS 2 , the total energy AE — AEj + AE t is min- 
imum, —3Nx 2 t 2 /8JS 2 . From this estimate, the system 
cants all the way to FM when the optimal 9 equals tt, 
which occurs at x = US 2 /\t\ ~ 0.025. 

Because the doubly-degenerate E g electrons hop 
anisotropically, a better canted solution exists. Let all 
spins with exp(iK ■ l) = 1 (K is (0,7r,7r)) tilt towards 
+x and all spins with exp(iK ■ I) = —1 tilt towards 
—x with angle 9/2, as shown in Fig. |l|. The magnetic 
energy cost is reduced to AEj — zNJS 2 (l — cos#)/6, 
but the delocalizaton energy lowering is also reduced 
to AE t = — 2n\t\ sin(0/2). This type of canting termi- 
nates in C-type antiferromagnetism when x = 2JS 2 /t, 
as shown in Fig. [p. Fig. || summarizes our results for the 
energies of uniform solutions versus x. 
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FIG. 2. Energy per electron of various doped states. Elec- 
tron concentration x is measured in unitsof 1//3 = JS 2 /\t\. 
Bold lines denote the most favorable uniform states. 
Non-uniform (polaron) states (discussed in Sec. VI; using 
F=0.25) are denoted by dot-dashed lines for the various val- 
ues of (3 indicated. Among different extended states, AFC lies 
lowest from f3x = 2 ~ 4, and FM takes over above /3x = 4 
(x ~ 0.025). For light doping, the inhomogeneous polaron 
states are stable. As x increases, uniform states take over. 
The transition from polaron-doped AFG to FM occurs at 
x ~ 0.045 (using T = 0.25, t = 0.75 eV, JS 2 = 4.74 meV, 
(3 ~ 158). Notice that the value of /3 determines the criti- 
cal doping x c of the phase transition, and even the transition 
type, as small (3 results in AFG to AFC transition, rather 
than to FM. But for small /3 the predicted x c can be large 
enough to exceed the light-doping approximation used here. 



IV. PURE SPIN POLARON EFFECT 

Experimentally, CaMnOa is insulating at small doping 
x. This rules out uniform solutions with electrons doped 
into delocalized states. The simplest picture is that each 
doped electron is localized on one Mn site, creating a lo- 
cal Mn 3+ with spin S = 2 in an AF background, for a 
spin excess of 1/2. This does not agree with the mea- 
sured saturated magnetization j^], which has been inter- 
preted in terms of local spin flipping (net excess spin of 
2) in the region 0.02 < x < 0.08. Local spin flipping 
leads to local FM regions (called spin polarons) around 
doping centers, allowing the doped electron to gain de- 
localization energy. In this section we discuss candidate 
localized ground states of a single E g electron, using the 
same Hamiltonian TL = TL t + ■ 

First we make a continuum (effective mass) approxi- 
mation in the spirit of Nagaev §]. Inside the local FM 
region, electrons hop like free electrons with inverse band 
mass 1/m* = (a/?i) 2 |i| given by Eqn.(||), with a the per- 
ovskite lattice constant, a ~ 3.73A This electron sits 
in a spherical well with infinite walls at radius R, whose 
depth is — 3|t|. G-type antiferromagnetism resumes for 
r > R. The ground state energy of the electron in this 
polarized region is then 

AE t = \t\(-3 + ^), (10) 

where the second term is the zero-point energy in the 
well. The magnetic energy cost from spin-canting is 

AE., = zJS 2 (AirR 3 /3a 3 ), (11) 

where the last factor is the number of Mn atoms in the 
cluster. Optimizing AE = AE t + AEj over R, the op- 
timal size cluster has a number of Mn atoms inside the 
sphere equal to 

4 /R\ 3 4 / ir\t\ \ 3/5 „„ 

This is close to the optimum 25-site symmetric cluster 
which we will find by exact diagonalization and show in 
Table I. However, we will also find a smaller asymmetric 
cluster of lower energy. 

It is interesting that continuum theory gives a lower 
energy if the region r < R is canted rather than fully FM. 
The A sublattice is fixed for all r, but the B sublattice is 
tilted toward the A sublattice by angle 9 for r < R. Then 
the optimum tilting angle is 180° (FM) for (3 = t/JS 2 > 
(3 C = (77r 2 /18)- 5 / 2 4z/7r ~ 220, while for smaller (3, the 
optimal tilting angle is sin(6*/2) = f3/f3 c . Our estimated 
value is [3 ~ 158, which gives (3/(3 c — 0.72. So the spins 
inside radius R are approximately 90° apart, and the 
optimum cluster increases to 31. 

Now we repeat the calculation using the true discrete 
Hamiltonian. First consider flipping only one spin. In 
this way, the spin-flipped Mn (the central site), along 
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with its 6 nearest neighbors, form a 7-site cluster with 
all 7 spins parallel. The cluster is invariant under trans- 
formations of the point group Oh . If the central Mn spin 
is kept unflipped, but instead the spins of those 6 nearest 
neighbors are flipped, a 25-site Oft-symmetric cluster is 
formed. Similar steps can be taken to obtain larger and 
larger symmetric clusters. We will not look at candidate 
states with canted spins except for one special case to be 
mentioned later. 

Each spin flip costs magnetic energy 6 x 2JS 2 . Since 
spins are parallel inside the cluster, the electron can hop 
among the E g orbitals of all the spin-aligned Mn ions, 
with a corresponding energy lowering from derealiza- 
tion. Table I shows numerical values of ground state en- 
ergy found by exact diagonalization of 2M x 2M Hamil- 
tonian matrices for symmetric clusters of M Mn atoms, 
ranging in size from M=l to 63. The ground states are 
all doubly degenerate (E g representation) because of the 
Oh symmetry of the cluster. 

Asymmetric clusters are also possible. If we flip one 
spin at site I and another one at site / + x + y, a 12-site 
cluster is formed. Starting from this 12-site cluster, there 
are four different ways to create a larger cluster with one 
more flipped spin, as shown in Fig. ||. Other possibil- 
ities are less closely packed, such as the 13-site cluster 
shown in the figure, and other 3-spin-flipped cases not 
shown here. The ground state energy of these examples 
are calculated, as shown in Table II. For our chosen val- 
ues of t and JS 2 , the most favorable spin polaron is an 
asymmetric 1 7-site cluster with three flipped spins. 

TABLE I. Ground state energy of symmetric clusters, with 
only spin polaron effects included. The last column uses 
t = -0.75 eV and JS 2 = 4.74 meV. 



cluster 



number of 
spins flipped 



energy gain from 
hopping (|t|) 



ground state 
energy (eV) 



1 

7 

25 
51 
57 
63 




1 

6 

13 
14 
19 





V3 = -1.732 
-2.330 
-2.449 
-2.380 
-2.600 





-1.242 
-1.406 
-1.097 
-0.989 
-0.869 




as3— 4 as3— 3 ^ 

FIG. 3. Different asymmetric clusters: as2 is the only one 
with two spins flipped; there are four types of as3 which have 
3 spins flipped. 



TABLE II. Ground state energy of asymmetric clusters, 
with pure spin polaron effect. 



cluster 



number of energy gain from ground state 



size spins flipped hopping (jtj 



energy (eV) 



12(as2-l) 
13(as2-2) 
16(as3-2) 
17(as3-l) 
17(as3-3) 
17(as3-4) 



-1.936 

-2 
-2.015 
-1.936 
-1.984 
-2.145 



-1.338 
-1.386 
-1.341 
-1.281 
-1.317 
-1.438 
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There are many possible variations with inhomoge- 
neously canted spins, of which we considered only one, 
a 25-site cluster with the 7 inner atoms canted rather 
than ferromagnetically aligned with the outer 18 atoms. 
This interpolates between the 7 and the 25 atom cluster. 
Without canting, the 7 and 25 atom cluster become equal 
in energy for the value \t\/JS 2 ~ 100, smaller than our 
preferred value of 158. It turned out that in the range 
94 < \t\/JS 2 < 104, the locally canted state was lower in 
energy than either the 7 or the 25 atom pure ferromag- 
netic cluster. 



V. LATTICE POLARON EFFECT 

The degenerate ground state, found in the previous 
section for symmetrical spin polarons, is Jahn- Teller (JT) 
unstable p4j. We now add to our Hamiltonian lattice 
distortions, controlled by the electron-phonon interac- 
tion. The only lattice degrees of freedom included are 
oxygen motions along the bonds to the nearest two Mn 
ions. If an oxygen moves along this bond a distance u, 
it gets closer to one Mn atom and farther from another. 
This will raise by gu (for the closer Mn) or lower by gu 
(for the farther Mn) the energy of any occupied Mn E g 
state of the type which points toward the oxygen (that 
is, the 3z 2 — r 2 state if oxygen motion in the z-direction 
is considered.) We use adiabatic approximation (oxygen 
mass — > oo ) and treat the oxygen distortions as classi- 
cal parameters. Each Mn ion is surrounded by six oxy- 
gens whose distortion amplitudes (us(H8/2), 6 = x,y, z) 
form basis vectors for a representation of Oh, namely, 
dig ®e. g ®t\ u . The vector representation t\ u contains the 
distortions u z (l + z/2) + u z (l — z/2), etc. The remaining 
3 degrees of freedom (Q z = u z (l + z/2) — u z (l — z/2), etc.) 
form basis vectors of a\ g (Q\ = (Q x + Q y + Q z )/V6) and 



H hl = -2s/2gJ2Qi(l) (4(1)^(1) + cl(l)c 3 (l)" 



(16) 



e g (Q 2 = (Q x - Q y )/2, Q 3 = (2Q Z -Q x - Q y )/Vl2), in 
Van Vleck notation [^5l . The oxygen distortions are lim- 
ited by a harmonic restoring force, e.g., —Ku z (l ± z/2). 

The lattice elastic energy and electron-phonon interac- 
tion terms of the Hamiltonian are |26| ] 

w L = — X)Mf+ U) 2 + u v(i + \vf + U S + \*)% 



(13) 



n cp = -^Y.^{i)c z (i)QS) 

+rotations to x, y directions]. (14) 
TLcp can be split into two parts, JT and "breathing" : 

'Hep = Hjt + Tibi (15) 
where, in Van Vleck notation, 



HjT = 2gJ2( 



4(0 



4(0 



Q 3 (0 
Q 2 (0 



0.2(1) 

-Q 3 (l) 



C2 (0 
c 3 (0 



The vector distortions (t\ u ) do not couple to E g electron 
states and therefore do not appear. 

A simple case shows how Ji jt splits energy degeneracy. 
Suppose a doped electron is localized at a single Mn site, 
with no hopping or spin flipping considered. The two 
E g states at that Mn are the only degrees of freedom 
for the electron and are originally degenerate. When the 
six surrounding oxygen distortions are considered, the 
degeneracy is lifted in Hjt: 



Hit - 2g( c \ c\)( 
= 2g ( at /Jt ) 



where 



Q 3 Q 2 
Q 2 -Qa 

-Q 
Q 



sm ■ 



C2 
C3 

a 



(17) 
(18) 



and (Qz, Qs) = (3(sin^>,cos0). The change in energy due 
to Hjt + H hl + H L is AE = -2V2gQi ± 2gQ + \2(Ql + 
Q 2 )/2. The energy splitting ±2gQ comes only from the 
JT term. The optimal distortions are Qi — \\[2gjK and 
Q = Ag/K, which give the maximum energy lowering of 
ground state AE = —6g 2 /K. 

The angle cf> does not enter AE. This continuous de- 
generacy can be lifted either by adding kinetic energy and 
quantization of lattice degrees of freedom (dynamic JT 
effect), or else by introducing higher order anharmonic 
terms in Tl cp [g7| pq |. However, later in the discussion 
of the 7-site cluster, we shall show how this continuous 
degeneracy is naturally removed by the increase of elec- 
tronic and lattice degrees of freedom. 

It is interesting to consider what would happen if spins 
were ferromagnetically ordered, so that magnetism does 
not assist localization. Then polaron formation can only 
occur through lattice distortion and is prohibited when 
the derealization energy per electron is larger than the 
polaron energy, i.e., when T = g 2 /K\t\ is less than a 
critical value close to 0.5. We believe that T is close 
to 0.25, so that in the hypothetical ferromagnetic case, 
polarons would not form. Antiferromagnetic confinement 
is needed before a lattice polaron effect occurs. LaMnC>3 
is different in this respect; its cooperative Jahn- Teller 
ground state makes polaron formation easier pi. 



VI. THE 7-SITE CLUSTER: SPIN AND LATTICE 
POLARON 

With the full Hamiltonian considered, the 7-site clus- 
ter turns out to be the most important one, as shown 
later in this section. Only a single spin is flipped. In- 
homogeneous canting will not be favored for this state. 
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By understanding its ground state algebraically, we learn 
several interesting aspects of the influence of H cp in the 
Hamiltonian. The dimensionless parameter T = g 2 /Kt 
characterizing the strength of electron-phonon coupling 
has a value near 0.25 for LaMnC>3, and we assume that 
the value for for CaMn0 3 is similar, i.e. 0.20 < T < 0.30 
|50| ^9|. We will measure energies in units of \t\, us- 
ing dimensionless lattice variables Q 1 = \j 'K/\t\ Q. The 
prime in Q' is suppressed from here on. 

To study the ground state of a single electron in the 
7-site cluster, a 14-dimensional Hilbert space is used, 
consisting of atomic E g orbitals from each of the 7 Mn 
atoms. Later when we turn on H op , we will find that for 
r small, fewer than 14 basis functions are needed for the 
ground state calculation. The 14-dimensional space can 
be decomposed into 7 irreducible representations of point 
group O h , namely, A lg @ A 2g ®T lu ®T 2u ®3E g . When lat- 
tice distortions are absent (T = 0,H cp = Hi = 0), these 
functions diagonalize the 14-dimcnsional Ht- The 3 E g - 
type basis functions (ijj 2 > ^3)1 (V*2) ^3) an d {^21 ^3) are 
degenerate separately with eigenvalues +V3|t|, 0, and 
— -\/3|i I . These 6 states, along with the A\ g state (with 
eigenvalue 0), will be the main states of interest when 
r 7^ 0. All other states stay absent as long as T is small. 

When r^O, H ep and Hi are turned on. Amplitudes 
of lattice distortions appear linearly in the matrix repre- 
sentation of H cp , and quadratically in that of Hi- Since 
these lattice distortions are of order W when T < 1, H ep 
and Hi can be treated as perturbations to Ht ■ In the 
following, all lattice distortion modes will be introduced. 
Then second order perturbation of the original ground 
states (ip 2 , ^3~) can be expressed in terms of these modes 
and shows that most of these modes do not participate 
in the ground state lattice distortion. This will in turn 
eliminate the need for considering electron states of sym- 
metries which couple to the absent modes only. Finally, 
from the form of the perturbed ground state energy, the 
pattern of its optimized lattice distortion will be derived 
and compared to the exact numerical result. 

There are 36 oxygens adjacent to one or more Mn 
atoms in the 7-atom cluster, so there are 36 distortion pa- 
rameters in H cp . These can be organized into sets of basis 
vectors for irreducible representations of the point group 
Oh, namely, 2>a\ g © a 2g © 5f i u © 3i2u © 4e g . The ai s -type 
modes are denoted as (qn, where i — 1,2, 3), and the e g - 
type modes are denoted as (qa,qa), where i = 1,2,3,4. 
The matrix elements of H op can be reexpressed in terms 
of these modes, and H L is simply £ roodes 9Ld cs / 2 - 

Second order perturbation theory shows that for small 
r, many of these modes appear only quadratically with 
positive coefficients and hence should be optimized to 
0. For larger T, these coefficients start to turn negative, 
and the corresponding distortions start to develop. The 
critical values are T = 0.443 for the appearance of the 
fi u and (743 distortion, T = 1.30 for t 2u , T = 0.819 for 
a 2g and q4 2 . It will thus happen that as T increases, 
the symmetry of the ground state electron wavefunction 



and lattice distortion pattern is gradually lowered by the 
successive appearance of these distortions. 

We therefore ignore the above modes for the actual 
range T w 0.25 ± 0.05. This eliminates the presence of 
electron states of A 2g , T\ u and T 2u symmetry. The re- 
maining modes are 3 sets of (qn, qa, qa), with i — 1,2, 3, 
of the inner, intermediate and outer oxygen layers, re- 
spectively. Degenerate first-order perturbation theory for 
(ij) 2 , 1P3 ) shows that the ground state energy of Ht splits 
into 



EV 



\t\ 
Qi 



1 



V5- Vf^QiT 
1 1 



(19) 



3«uH 


- gfti H 


- 12931, 




1 
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Thus at first order, there is still a continuous degeneracy, 
in the sense that the energy JT splitting depends only on 
Q = y/Q'i + Q'i, not on <f> = t&n- 1 (Q 3 /Q 2 ). 

To include 2nd order perturbations, for simplicity, we 
treat analytically only the distortion modes from the in- 
ner oxygen layer, (911,912,913)- The number of related 
electronic states is now reduced to 5, namely, they are 
{Ai g , , ip£) and (ip 2 , V';}")- I n this 5-dimensional sub- 
space, the Schrodinger equation to be solved can be ex- 
pressed as follows (not including Hi,, which is always 
proportional to the identity matrix): 
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where * J ee (tJj 2 , ip 3 ), * 
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an effective Hamiltonian H e s = H{ + H„ n — h^H 11 
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(A lg , ip£ , ^3 )■ We obtain 
t + ^cp - hl(H t + 
... . ... ■ Since E should be very close to 

— V3|f| for small perturbation, we take E as — V3|f| x 
1(3x3) and ignore H^ p in the denominator. The H e s 
obtained this way shows that the degeneracy of (tp 2 , %p 3 ) 
is now lifted to become (including Hi) 

4 4 32 a2 2 8a/6 , 2 2 
e = 3P + Y A p + ~ Ap + p 

+ ^p 3 (4V2A + V3) cos (36>), (21) 
where we introduce the notation 



B = —912 = psmO, C = —9i3 = pcos( 



6 



We see that cos 36* should be +1 to minimize the ground 
state energy. Hence the degeneracy of ground states is 
not continuous anymore and has become 3-fold. This 
feature agrees with the numerical result. The ground 
state electronic wavefunction of 6 = is shown in Fig. 
^. A rotation which brings z to x (y) will generate the 
wavefunctions of 6 = 2ir/3 (9 = 4tt/3). 




FIG. 4. One of the possible ground state wavefunctions in 
the 7-site cluster. The corresponding lattice distortion pat- 
tern has qi and q3 components. 



6 oxygens surrounding the central Mn. Displacements of 
oxygens in the intermediate and outer layers are at least 
10 times smaller. As predicted perturbatively, t\ u , t 2u , 
o-2g, <?42 and (743 distortions are absent. 




VII. NUMERICAL RESULTS 



FIG. 5. Numerical results for ground state energy(_E 9S ) of 
symmetric clusters 



Now we consider the full Hamiltonian H = Ti t + Tij + 
7i C p + "Hl for all clusters examined in Sec. IV. FM spin 
alignment facilitates hopping and causes energy lower- 
ing from derealization, therefore encourages the polaron 
to grow. However, as T increases, the JT splitting will 
become the dominating influence on ground state en- 
ergy. Greater localization enhances the JT energy lower- 
ing, which increases localization of the electron. On the 
other hand, Tij and 7Yl serve as penalties for spin mis- 
alignment and lattice distortions. Our numerical studies 
find the optimal resolution of the competition between 
these effects. The 7-site cluster becomes favored from 
r > 0.18, as shown in Fig. ||. and Fig. ||. For all clusters 
considered, when T becomes large enough, the size of the 
ground state wavefunction shrinks to become that of the 
7-site polaron, as shown in Fig. [?|. Enlarging the size of 
the FM cluster to enhance derealization energy is then 
disfavored by the strong lattice polaron effect. It is also 
clear that for extremely large V, the 1-site lattice polaron 
will be the only form of electron state which exists. 

Our numerical calculation predicts the lattice distor- 
tion pattern of the 7-site ground state. For T = 0.25, 
t = 0.75 eV, and K = 27.2 eV/A 2 (obtained from Ra- 
man scattering in LaMnOs), the ground state shown in 
Fig. | has oxygen displacement parameters qu ~ 1.1, 
<7i3 ~ 0.78 and = 0, which gives an outward displace- 
ment of 0.38A in the z direction, and smaller outward dis- 
placements of 0.036A in both x and y directions, for the 
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FIG. 6. Numerical results for ground state energy (E gs ) of 
asymmetric clusters 
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FIG. 7. Numerical results for ground state radius of sym- 
metric clusters. The radius is measured in unit of lattice 
constant a(the Mn-Mn separation distance), and is defined to 
be [y\(rj— < f >) 2 ?/)f] 1,/2 , where the index i runs over all 
Mn sites inside the cluster, and < r >= y"V . fiip^. 

VIII. DISCUSSION 

When doping x is non-zero, one should ask whether 
localized polaron solutions will distribute homogeneously 
or will attract, causing phase separation [|T). We have 
not addressed this issue, which requires a more compli- 
cated calculation with additional Coulomb parameters in 
the Hamiltonian. Experiment is consistent with a con- 
centration interval up to x = 0.08 where polarons are 
homogeneously distributed. Our model shows that by 
this concentration, randomly-distributed 7-site polarons 
will overlap significantly. However, our model also shows 
that at concentrations x > 0.045, polarons should be un- 
stable relative to an undistorted ground state with FM 
spin order and metallic conduction by the doped elec- 
trons. This effect is not seen in experiments. Apparently 
alternate ground states, possibly involving organization 
of polaronic distortions, occur and enable the system to 
remain non-metallic. 

Without additional physics (such as defects), our 
model cannot account for the observation |j| that La con- 
centrations x less than 0.02 yield less excess moment than 
expected from 7-site polarons. 

Our model describes the competing spin and lattice 
polaron effects. We believe that it contains the main 
features needed to describe the system. A test would 
be measurement of the oxygen displacements which our 
model predicts. The model omits non-adiabatic phonon 
effects, spin quantization, temperature, and polaron- 
polaron interactions [p2[. For higher doping levels or 



T > 0, these may have a larger influence and present 
challenges which could be worth pursuing if experimen- 
tal guidance improves. 
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